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Abstract 

Two puzzling facts of our time are the observed patterns in the fermion masses and mixings 
and the existence of non-baryonic dark matter, which are both often associated with extensions 
of the Standard Model at higher energy scales. In this paper, we consider a solution to these 
two problems with the flavour symmetry A 4 x Z 2 x Z 2 , in a model which has been shown before 
to explain large leptonic mixings with a specihc texture. The model contains 3 generations 
of S'17(2)i-doublet scalar fields, arranged as an A 4 -triplet, that spontaneously break the elec- 
troweak symmetry, and a “dark sector” of Z 2 -odd fields, containing one Majorana neutrino and 
an A 4 -triplet S'17(2)4,-doublet scalar field, the lightest of which provides a candidate for dark 
matter. 

Concerning the Z 2 -even scalar fields, compared to the Standard Model, we predict additional 
fields with masses at the electroweak scale. We therefore investigate present phenomenological 
constraints from lepton flavour violation experiments, obtaining a lower bound on the extra 
scalar masses of 140 GeV. Furthermore we consider the oblique parameters, Higgs boson decay 
properties and possible flavour violating signals at the LHC. 

Concerning the “dark sector”, we study bounds from dark matter search experiments and 
identify the parameter space of the dark matter candidate that is compatible with the observed 
relic density. We hnd two allowed mass ranges for the dark matter within which the experimental 
constraints can be accommodated: the low-mass range is from 47 GeV to 74 GeV and the high- 
mass range is from 600 GeV and 3.6 TeV. 


1 Introduction 

The Standard Model (SM) continues to be an extremely successful description of particle physics. 
The scalar resonance that was observed by the ATLAS and CMS collaborations [^[^, with mass 
125.7±0.4 GeV is consistent with the SM Higgs boson. Nevertheless, there is still room for 
extensions of the scalar sector and despite the discovery of a SM-like Higgs, the motivation for 
physics beyond the SM remains with the many unresolved puzzles that are not addressed by the 
SM. One of these puzzling facts is the existence of non-baryonic dark matter (DM), which is about 
five times more abundant than baryonic matter [^. Up to now it is undetected by direct detection 
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experiments and potential hints for DM signals in indirect detection experiments are not always 
conclusive due to astrophysical backgrounds, see e.g. refs. [5||^. A new unbroken symmetry seems 
to be of fundamental importance to explain the DM properties; it prevents the DM particle from 
decaying into SM particles, and thus stabilises it while also explaining the suppressed interactions 
with the SM particles. A simplistic “bottom-up” approach to the DM problem is given by the class 
of one-particle extensions of the SM. A more “top-down” approach often generates an entire “dark 
sector”, as is the case e.g. for supersymmetry, extra dimensions and grand unified theories. 

Another puzzling fact that remains unresolved within the SM is the observation of regularities 
in the fermion sector of the SM. The mere existence of three fermionic generations cannot be 
addressed or understood within the SM let alone the particular pattern of the fermion masses 
or their mixings. A very peculiar feature is also found in the lepton sector, which shows at the 
same time tiny neutrino masses and large mixings between the generations. This is a pattern that 
is different from what is observed for the quark sector. 

An extension of the SM by one or more additional Higgs doublets (see e.g. ref. 11 for a 


review with the focus on two-Higgs-doublet models) can address both of the above mentioned short¬ 
comings of the SM. One particularly intriguing possibility is to introduce the additional doublets 
as a multiplet of a flavour symmetry such as A 4 [12-15 , see also ref. 16 and references therein. In 


such models with a flavour-symmetric extension of the scalar sector the Higgs phenomenology leads 
to important direct constraints, while indirect constraints are given by fermion flavour violating 


processes or simply due to the requirement that the fermion masses and mixing are viable 17 -19 


In this type of A 4 model it is possible to account for an excess in /rr events that was found by the 
CMS collaboration |20l|21j, which is interpreted as lepton-flavour-violating decays of the SM-like 


Higgs boson, see for instance ref. 22 -24 


In this context, a very interesting “top-down” approach is given by the connection of the DM 
with the flavour sector where the DM candidate can be stabilised by the same symmetries that are 
introduced to explain patterns of leptonic mass and mixing, see e.g. refs. (2^29j for examples with 
A 4 symmetries. In flavour models with multiple scalar fields that have masses on the electroweak 
scale, i.e. O(TeV), there is the promise of rich collider phenomenology that can potentially be 
explored with the next LHC run, at 13 TeV. In these models the DM candidate is typically part 
of a larger “dark sector”. 

In this paper we study a multi-Higgs and multi-component scalar DM model that entails a 
rich phenomenology, due in particular to the additional scalars having masses on the electroweak 
scale. We show how existing measurements constrain the model and how it can be further tested 
in ongoing and future experiments, such as the next run of the LHC and dark matter search 
experiments. The model is inspired by the A 4 model in ref. [16] , where the neutrino masses are 
generated radiatively and the observed large leptonic mixings can be reproduced. 

This paper is structured in the following way: we review the relevant features of the model under 
consideration in section in particular the scalar potential and the leptonic Yukawa couplings. 
We investigate phenomenological constraints on the Z 2 -even scalar fields in section in particular 
from lepton flavour violating effects and Higgs decay properties. In section we consider the 
phenomenology of the Z 2 -odd scalar fields, where we in include constraints from astrophysics and 
direct searches at colliders. We conclude in section In appendices we include the details of A 4 
and also the A4 x Z 2 x Z^-invariant scalar potential used in this paper. 


'*An argument has been presented in [^ which connects the SM gauge group with the number of fermion genera¬ 
tions. 
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Table 1: Field content with charge assignment of the model. Note that x is real. 


2 The Model 


The model considered in this paper consists in a multi-Higgs extension of the SM, with flavour 
symmetry A 4 ( 8 ) Z 2 ( 8 > Z 2 [16] . The field content is summarised for convenience in table with its 
symmetry assignments. 

We embed the left-handed SM leptons into an A 4 -triplet and add one singlet Majorana 

neutrino Nr that is odd under Z 2 . The scalar sector contains two S? 7 ( 2 ) 2 ,-doublet A 4 -triplet fields 
(1 2 ) 

, one of which also carries a Z 2 -charge, and another A 4 -triplet field X; which is charged under 
the symmetry Z^. The Yukawa part of the model Lagrangian for the lepton sector takes the form 


C-Y = Ve ^eR + Ufj. fJ-R + Vt _^JR 

+ Vu ^ Nr + MnNrN^r + h.c., 


( 1 ) 


with = icr 2 {k = 1, 2). The subscripts 1 , 1 ', l" denote the projection of the correspond¬ 

ing A 4 singlet in the product of the two triplets. For a brief summary of the A 4 product rules, see 
Appendix To be explicit we define the following decomposition of the scalar fields: 








-|- hj 


i hj 


with j = 1 , 2 ,3 . ( 2 ) 


From the Yukawa terms in eq. Q the charged-lepton Yukawa matrices in the flavour eigenbasis 
^= 1 , 2,3 for each component of follow: 


/ Ue Vu yT\ /O 0 0\ /OO 0\ 

Fi = [ 0 0 0 I , Y 2 = f 2 /e wyJ , Y 3 = f 0 0 0 ) 

\ 000 / \0 0 0 / \ye voVfj, 


(3) 


Due to the Z 2 symmetry, the field Nr couples exclusively to the A 4 -triplet iS'17(2)-doublet 
The Yukawa couplings for the neutrinos involving the components and Irj are given by yy 6ij. 

When the A 4 -triplet SU (2)^.-doublet develops a non-zero vacuum expectation value (VEV), 

the Yukawa matrices in eq. ^ generate a mass matrix for the charged leptons. Without breaking 
the Z 2 symmetry, the masses for the left-handed neutrinos can be generated radiatively, via 
and (x). The corresponding mass matrix of the neutrinos is sensitive to both VEVs, in contrast 
to the mass matrix of the charged leptons. As was shown in ref. |16|, in this case viable leptonic 
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mixing is possible if the two VEVs {x) have a specific alignment in A 4 -direction and the 

magnitudes of the VEVs exhibit a hierarchy 




=: r <C 1 . 


(4) 


The most general scalar potential which is invariant under the SU{2)l x A4 x Z2 x Z 2 -symmetry 
can be expressed as: 


V = $(1), $(2)) + v^^ix) . (5) 


The three sub-potentials on the right-hand side of eq. ([^ are shown in full in Appendix We 
find that for a suitable choice of the parameters, the sub-potential can indeed cause the real 
scalar field x to develop the VEV (x) and break the A 4 symmetry. The sub-potential consists 
in the interaction terms between the x fields and the = 1,2 fields. It also adds terms 

proportional to (x) to the mass matrices of the fields i = 1,2. The sub-potential allows 
for electroweak symmetry breaking and adds interactions between the scalar fields and 
This sub-potential can lift the mass-degeneracy between the components of the three electroweak 
doublet generations contained in 

The special directions in A 4 -space of the VEVs (x) and are referred to as VEV alignment. 

The particular VEV alignment that was used in in order to generate large leptonic mixing is 


= 0 , 



(X) = ^(l 0 -1), 


with uew = 246.22 GeV <C Vy. being the electroweak vacuum expectation value, cf. eq. (1]). The 
horizontal direction is the A 4 direction where the vertical direction (displayed for shows the 

SU{2)l components. Note that the field does not develop a VEV, such that the Z 2 symmetry 
remains unbroken. For the potential in Appendix [B| the direction of (x) is a stationary point but 
it is not a true minimum. We discuss this issue in section]^ and in Appendix [C} 

With the VEV ($( 2 )^ as in eq. ®, the mass and mixing matrix, respectively, for the charged 
leptons follow from eq. (§: 
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such that is diagonal and positive. In turn, the structure of the light neutrino mass 

and mixing matrix, respectively, is dominated by (x) and takes the form 
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where a and b are specific loop functions that depend on the Yukawa couplings, mass of the RH 
neutrino and the masses of the Z 2 -odd scalars, and P^, is a diagonal phase matrix. V,^ and 
combine to a viable PMNS matrix as described in 
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We remark that in this model the assignment of quarks under A 4 was not considered explicitly. 
In A 4 -symmetric models where the scalar S'[/(2)i-doublets transform as triplets, the quark masses 
and mixings can be generated by embedding the three generations of up and down quarks into 
A 4 -triplet fields. Thereby either the left-handed S'17(2)i-doublet Ql, the right-handed SU{2)l- 
singlets ur and dn, or both, can be defined to transform as A 4 -triplets. This type of assignment 
with quark generations embedded in triplets is usually considered in grand unified models featuring 
A 4 to explain leptonic mixin g [30] ( see |31| for a recent example), or otherwise to obtain relations 
between quarks and leptons (T2l|33| . 


3 Z 2 -even Scalar Sector 

3.1 Physical Field Content 

We observe that the VEV of the Xi fields breaks the Z 2 symmetry, such that they mix with the 
fields. The mass matrix of the full scalar sector is given by 


M(Xkj = 


d^V 

dXidXi 


X = h,h,h^,rj,fi,r]^,x 


(9) 


with the helds h, rj from the dehnition in eq. Q and where we suppressed the A 4 indices. The 
charged helds 2 3 liiiear combinations of the charged ^^= 123 - physical scalar helds 
Hi, i = 1,...,6 and llj=i, 2,3 are linear combinations of the havour eigenstates hk,hk and Xk, with 
k = 1,2,3. They are obtained from diagonalising the mass matrix in eq. Q, which is block 
decoupled from the helds ??i=i, 2,3 due to the unbroken Z 2 symmetry. In the limit r = 0 the mixing 
between the h, the h and the y vanishes and the transformation between havour basis and mass 
basis can be expressed as 


H = TZhX 




X = {hi,h2,h3,xi,X2,X3y 

x= {hi,h2,h3'^ 


( 10 ) 




X={hf,hf,hff , 

with the helds h, x from the dehnition in eq. the 6x6 matrix TZh being given by 


TZh = 


7^<j> 


( 2 ) 


( 11 ) 


and the rotation matrices ^^d H 




0 - 


1 \ 
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( 12 ) 
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The mass matrix defined in eq. Q contains some eigenvalues that vanish in the strict limit r = 0, 
therefore we present the eigenvalues to leading order in r: 

'^Hi = 4 (-^ 1 ,( 2 ) + ^ 2 ,( 2 ) + A 3 ,( 2 ) - lai+k) ^EW > 

~ ^ (^4,(2) + A5,(2)) "ylw ) 

= 2 (72,(2) ’ 

rriH^ = 4(3^2 - 4(ii)u^ , 

= 2 (4di - 3(^2) v‘^ , 
rriH^ = 4 (4di + <^ 2 ) v'^ . 

For the charged fields we get the following eigenvalues, using the same rotation matrix TZ^ 

= 0, 


(13) 


( 2 )’ 


m 


— (A5,(2) “ 2(A2,(2) + A3,(2))) ) 


fj± 2 (72,(2) I'x • 


(14) 


Finally we hnd for the eigenvalues for the fields Ffi=i, 2 , 3 , again using TZ^ 

A 
A 

,2 


(2)- 


rrif, = 0. 

Hi 


"^^2 — 2 (A5,(2) - 2A2,(2)) VeW 


(15) 


m 


Hs ^ ^2,(2) ^x 


Among those, the neutral field Hi and the charged field with mass eigenvalues zero can be 
identified as the Goldstone bosons. The held Hi is the SM-like Higgs boson of our model, that we 
can identify with the scalar resonance at 125.7 ± 0.4 GeV [^. It couples to the mass eigenstates 
of the leptons via the Yukawa matrices 



1 - 


^1.(2) 


4di + d2 


(16) 


which is proportional to the SM Higgs boson couplings with a reduction factor proportional to the 
small parameter r. 

The second generation of Z 2 -even helds, namely H 2 , H 2 , have masses proportional to the 
VEV of and thus also reside at the electroweak scale. They couple purely off-diagonally to 
the leptons, via the Yukawa matrix 


1 


0 y^uj‘^ VrUJ \ 

H 2 = A I ygcj 0 yrU}'^ , 

^2/eW^ VliUJ 0 / 


(17) 


with oj = The Yukawa matrices for the second generation of the Z 2 -even scalars are given 

by 

= i , {Yf),. = V2 , (18) 
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(a) T ^/ie (b) /i —?► 67 (c) e“''e —> r+/j, 

Figure 1: Lepton flavour violating processes, mediated by the fields H 2 and H 2 . 


while the analogous Yukawa couplings for the third generation can be obtained using the substitu¬ 
tions U! —)• —iw and —)• iaj^ in eq. © and eq. ( [I8| ). The fields iL 4 , 5,6 do not couple to fermions 
at all, as they originate from y. Note that beyond leading order Hi mixes only with one state 
contained in y (namely Hq) leading to the shown reduction factor compared to the SM and to the 
fact that all vanishing entries in Yukawa matrices also vanish for r 7 ^ 0. 

Notice that the squared masses of the iL 4 and in eq. (13) are not independent for r = 0: 

This is due to the VEV (y) being a saddle point solution rather than a minimum 
of the sub-potential V^-^. One can impose the relation di = 3 / 4^2 to avoid a tachyonic mass. This 
restrictive region of parameter space then allows (y) to become a flat minimum where, however, 
even for non-zero r one additional zero mass eigenvalue is present. We present a simple solution 
to this issue in Appendix [C| which does not affect the discussion of the phenomenology. For this 
reason, in what follows we shall consider the mass-scales of the iLi >3 fields to be 0(v^). 


3.2 Phenomenology 

The Z 2 -even scalar fields of the considered model exhibit the peculiar feature of having two gen¬ 
erations of fields that couple completely off-diagonally to the charged leptons. In this section we 
investigate the subsequent constraints from experimental bounds on lepton flavour violation, from 
electroweak precision data and from the decay properties of the 125.7 GeV SM-like Higgs boson. 


3.2.1 Lepton Flavour Violation 


Lepton flavour violation (LFV) occurs at tree level in our model, due to the Yukawa couplings of 
the second generation Z 2 -even scalar fields {H 2 , H 2 and H^) - the third generation in principle 
contributes similarly, but is assumed to decouple due to its mass of 0{v^) in the following. It is 
interesting to note that the typically very constraining process fj. —>■ eee is not mediated at tree 
level by H 2 or H 2 , since both fields have no flavour-conserving couplings to the leptons. Instead, 
the following processes are relevant: 


e 


/i , fj. e'^e , and e . (19) 

For these processes we depict an example diagram in fig. (a), with the flavour-violating fields H 2 
and H 2 . The leading contributions to the lepton flavour violating tau decays of eq. (19) are given 
by 
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Figure 2: BR(/x —)• 67 ) as mediated by the fields H 2 and H 2 - The field contributes negligibly. 
The dashed line corresponds to the experimental constraint BR(;U —)• ey) < 5.7 • 10“^^ li^. 
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( 21 ) 

( 22 ) 

(23) 


where we defined Rr^x = BR(r —>■ X) /BR(t —>■ e Vei^r) and the lepton masses = 0.105 GeV 


and rriT- = 1.77 GeV. All formulae were obtained using FeynArts [34] and FormCalc [^,36 . The 
process r“ —)• ii~e^e~ with the experimental upper bound being 0(10“®) results in the very mild 
constraint on > 12 GeV. 

Recent and very stringent constraints on rare charged lepton flavour violating decays come from 


aoration 
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the MEG colla 
depicted in hg. 
can be approximated as (for not too large A): 


The contribution from the helds H 2 and H 2 to those rare decays is 


^ (b). With the definition + A/2, — A/2, the contribution 


BR(;U —)> ey) = 


am" 


AM* - eOAM^ + 227A2 


3847rM8 L 


+ 24Alog(^) (-2M2 + 




15A + 6Alog(^))j 


(24) 


This was obtained using the formulae in ref. 38 


The charged field also contributes to the process, but according to the corresponding 


Yukawa matrix in eq. (18) its contributions to eq. (24) are suppressed by the tiny ratios me/m 
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Figure 3: Cross section of the dominating pro¬ 
cess contributing to the production of H 2 pairs 
at the LHC, for ^/s = 7 TeV. The dashed black 
line denotes the particular model parameter set¬ 
ting of m^± = mH 2 ) while the solid black line de¬ 
notes the irreducible contribution from the gauge 
interactions of the H 2 field. 


and m^/m-r- For the case ~ '^H 2 ' experimental bounds on the branching ratio of BR(/r —?■ 
67 ) < 5.7 X 10“^^ then requires mH 2 > 143 GeV. We note that when the model saturates the MEG 
bound, i.e. produces a branching ratio of // —>• eq comparable with the current sensitivity limit, the 
branching ratios for the similar processes in r decays are given by 

BR(r —> ^ 7 ) = 2.1 X 10“^® and BR(r — )■ 67 ) = 1.5 x 10“^® . (25) 


Bounds on LFV have been established at LEP-II, for instance by the OPAL collaboration 
Eor centre of mass energies from 189 to 209 GeV the following limits were obtained: 
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a{e^e —>• l(!) < 166 fb , 


(26) 


for W = fiT. The limits on e/r are slightly more constraining, however such final states cannot be 
produced at leading order within our model. The leading diagram involving H 2 , H 2 which leads to 
LEV probed at LEP-II is depicted in fig.j^(c). Eor = 140 GeV, the contributing scattering 

cross section is of order yoctobarn (i.e. 10 “® femtobarn, also referred to as a shed), which is far 
below the experimental sensitivity. 

It is worthwhile to consider LEV signatures of the Z 2 -even scalars at the LHG. The fields H 2 
and H 2 can be produced via vector boson fusioij^and decay into lepton flavor violating dibosons: 

PP 3 3 H 2 H 2 , with LEV decays H 2 ^ , (27) 


where j is an energetic jet, and the same process exists for the replacement H 2 —>• 772 - The 
dominant sub-process is given by ud —)■ duH 2 H 2 (and analogously for H 2 )- We show its cross 
section in hg. for ^/s = 7 TeV. The numerical values were obtained with GalcHEP 


GTEQ 6 parton distribution functions 41 
£ = 20.3 fb-i 


40 and the 


With the present integrated luminosity at the LHG of 
a total number of LEV events of ~ 10 can be expected for each, the H 2 and the 772 - 


The combinations of lepton flavour thus produced are er and //r. Note that the expected number 
of events for the lepton-flavour combination e/i is suppressed by the ratio 

Recently the direct search for lepton-flavour-violating Higgs boson decays of the GMS collab¬ 
oration yielded an excess of 2.5a in the Higgs-to-^u-r channel [^,21 . The LEV excess resides in 
the energy bin from 120 to 140 GeV and was fitted as an LEV decay of the SM Higgs boson. It is 


^Note that the interactions of the scalar sector with the quarks has not been defined within our model. 
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straightforward to consider this LFV excess being caused by the H 2 itself, which implies 
to be at the lower bound of 140 GeV, which implies in turn that MEG is about to observe the first 
// —)• ey events. Note that in order to account for the observed /i —)• /ur excess, more than 0(10) 
LFV events are necessary - the 10 events estimate presented above was considering only vector bo¬ 
son fusion processes. Interestingly, the 2-jet LFV events reported in are 4.1 ± l.l|^r^ ± 0.9|^T-e, 
which is around the order of the expected number of vector boson fusion events in our model. 


3.2.2 Effects of a Perturbation of the VEV Alignment 


The bounds above apply strictly speaking only if (73 ( 2 ) = < 74 ,( 2 ) = 0 in (see eq. (55)). While 
this may be possible at the classical level, quantum corrections, e.g. the ones induced by the terms 
in Appendix]^ can render these parameters non-zero, which deforms the VEV of as shown 
in eq. , to 


/^( 2 )\ 0 0 0 
\ / Ve VI1 + 25 l-5 + i\/35' 


(28) 


In the following, we assume the parameters 6 and 5' to be small in order to allow for a perturbative 
treatment. They are given to leading order in r and ( 73 ,( 2 )) *^ 4 ,( 2 ) by 


^ ^ '^3,(2) + <^4,(2) 1 

2\/2 (A4,(2) + ^5,(2)) ^ 

S' ~ ^3,(2) ^5,(2) + <^4,(2) (^4,(2) + ^ A5,(2)) 
2-v/2cJ2,(2) (A4,(2) + A5,(2)) 


(29) 

(30) 


where we assumed < 73 ,( 2 ) and < 74 ,( 2 ) to be real for simplicity. The requirement of 5, 5' being small 
demands the parameters < 73 ,( 2 ), < 74 ,( 2 ) to compensate a factor 1/r. Furthermore, due to r <C 1, the 
subleading perturbation proportional to 5' can be neglected. 

With the perturbed VEV alignment in eq. (28), the charged lepton mass matrix of eq. 0 is 
not diagonal after rotation with and requires additional rotations 9ij ~ 0(6) in the left-handed 
sectoij^ and 9ij ~ 0{6)yi/yj in the right-handed sector. As a result the zeros in the Yukawa 
matrices Yi, Y 2 and Y 2 (in the mass eigenbasis) are corrected to (up to leading order in 5) 


/ ye 0(6) yf, 

Y^ = —\0(6)ye y, 

\0(6)ye 0(6) y^ 

and 

1 (0(6) ye 

V yeW 

with analogous structure for l 2 - 


0(6) yr\ 
0(6) yr\ 

Vt ) 


0-l,(2) 

" \ 

(31) 

Ml + d2 

T 

0(6) y^ 

yrUJ ^ 

0(6) yr; 



(32) 


®While this might change the conclusions made for the flavour model this paper is based on, we will not go into 
more detail on this particular issue and assume it non-problematic. 
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The perturbed Yukawa matrices above allow for additional flavour-changing vertices for Hi 
and flavour-conserving vertices for H 2 and H 2 , such that non-zero contributions to the processes 
(.i —)• arise. These can be estimated as 


BR(£i ^ Uj) 


2 2 
— _ I 

4 ^ ' 

rrii 


BR{ii IjvU) ' m% 

where H is the scalar boson mediating the decay as described by diagrams of the form of fig. 
(a). Assuming H = Hi with ~ 126 GeV, the most stringent bound on 6 comes from r —)• 3/i 
and is given by 5 < 30. Since J < 1 in order to allow for the perturbative treatment, these kind of 
processes do not give meaningful bounds on the misalignment at present. 

Likewise, when considering the branching ratios BR(.^j —)■ Ij'^) as described by fig. (b) and 
similar, we find further contributing diagrams. Compared to the unperturbed VEV alignment, 
only the two processes r —)■ //y and r —)■ ey receive potentially signihcant contributions, as it is 
now possible to have tau leptons in the loop, with large Yukawa couplings to the scalars. Said 
contributions to the branching ratio are “enhanced” by a factor of m^/m^ 0((5^) ~ 280 5^. Com¬ 
paring the result for an unperturbed VEV from eq. (25) with the experimental bounds of 0(10“®), 
it shows that no bounds on 5 can be obtained from this type of processes either. 

Physical processes where the perturbation of the VEV alignment could lead to observable 
effects, are the LEV decays of the Hi boson. In particular the branching ratio BR(/i —)• /rr) can be 
estimated as 0(5^)BR(/i —>■ rr). This allows to use the observed excess of h —)■ /rr events at the 
LHC with a branching ratio of 0.84i’iQ3y% to determine 6 ~ 0(0.1). In this case, it is not necessary 
any more to have the masses of H 2 and H 2 close to the MEG bound. Eor a more comprehensive 
analysis of the effect of a perturbed VEV alignment in a similar framework, also considering the 
h ^ jjLT excess, see ref. (^. 


3.2.3 Oblique Parameters 


The difference between the masses of the real, pseudo-real and charged fields from the previous 
section can be constrained via their contributions to the oblique parameters S, T, U. This has 
been studied for instance in refs. 17,42 . The strongest constraint comes from the T parameter. 
A recent fit of the SM parameters to a global data set yields (4^ 


r = 0.09 ±0.13. 


(34) 


In order to assess the contribution to the T parameter we consider a generation of a neutral scalar 
field S, a pseudoscalar field S and a charged scalar field S^. We define the masses of the three 
fields by 

ms± = pms and rrig = prug ■ (35) 

We can approximate the contribution of one generation of scalar helds, S, S, S^, to the T- 
parameter by 

^ [{p - I? -{p- 1)' - (P - pf] ^ , (36) 

G-TTsfy 

with siY being the sine of the weak mixing angle and mw the mass of the W boson. We observe 
that for p = p OT p = 1 the contributions of the three fields cancel exactly. Furthermore, it is 
interesting to note that STg > 0 for p > p and 5Ts < 0 for p < p, which allows the possibility of 
cancellations between the contributions from different generations of scalars. 
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Figure 4: Higgs to diphoton decay width, normalised over the SM prediction, in the presence of 
one extra charged scalar boson of mass mg± and coupling to the Higgs boson The green 

line denotes the combined experimental precision from CMS |44|-l46| and ATLAS 47,48 . The 


blue line denotes the maximum deviation for |u; 5 ±| = 1 , the red line indicates the limit set by the 
LUX collaboration on the spin-independent-DM-nucleon cross section, which may be related to the 
Higgs coupling of the charged “sibling” of the DM particle, see section 4.3.3 


3.2.4 Diphoton Decays of the Higgs Boson 

The theory prediction for the process Hi 77 , where Hi is the SM-like Higgs boson in our model, 
is given by: 


r — 


128\/27r3 


Aw{tw) + E NcjQjApiTf) + ^ QxAxirx) 


fermions 




(37) 


with Tx = m‘jj^/4:mx, Nc^x being the number of colour states (3 for quarks, 1 for leptons), Qx is 
the electric charge of the particle in the loop, and the amplitudes A{t) are functions that depend 
on the spin and couplings to the Higgs of the particle running in the loop. 

The amplitude Ax denotes contributions from the electrically charged non-SM scalar fields, 
in our model given by X = with electric charge \Qx\ = 1- We define the fields 

Af=i 2 3 ™ section 0 below, as the mass eigenstates of the Z 2 -odd helds rj^Li 2 3 - The amplitude Aj 
for the field / = F,W, S being a fermion, vector boson or scalar, respectively, can be expressed 
as 


with the function 


Af{t) 

= (r + (r - l)/(r)) 


(38) 

Aw{t) 

= -;^(2T2 + 3T + 3(2r-l)/(r)) , 

(39) 

As{t) 

= - fir)) ] 


(40) 

fi'r) = 

I arcsin^\/r 

r < 1 

2 

r > 1 

(41) 
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Note that this definition of the Af includes the SM relation between the mass of a particle and its 
coupling to the 125.7 GeV Higgs boson Hi. 

We start by considering a generic electrically charged scalar boson S^, with mass mg± and an 
interaction term with the SM-like Higgs boson of the form; u!g± vew Hi . The contribution 

of the scalar field to the SM-Higgs-to-diphoton decay width can be expressed as 


, , 246 GeV , , 

-^s± G-s±) = wgi- As{ts±) , 

ms± 


(42) 


where the product on the right-hand side of eq. (42) takes into account the deviation from the SM 
relation between coupling and mass of a particle. 

We show the ratio of the decay rate T.^.^ - for the case where the SM is extended with one 
extra charged scalar field - over the SM prediction T®^ in fig. 4l Notice that an enhancement 


77 _ 

requires positive interference with the top quark, which occurs for ujg± <0. It shows that even for 
the coupling between the scalar and the Higgs boson being order unity, the deviations from the 
SM prediction are compatible with the experimental result for mg± > 200 GeV. For mg± on the 
TeV scale and with couplings order one the induced deviations remain on the percent level. This 
is beyond the reach of the LHC even with the luminosity upgrade. 

As we have mentioned above, in our model the electrically charged scalar fields relevant for 


and 2,3 


(see section for 


the diphoton decay width of the Higgs boson are the fields 77)^2 
the definition). In fig. we include the LUX constraint on the modulus of the coupling between 
the (electrically neutral) field and the Higgs boson I^dmI) which might be related to uf, the 
coupling between Hi and Af. 

The deviations from the SM prediction due to charged scalar fields with masses around the 
TeV scale could be studied for instance at future lepton colliders. The following precision can be 
achieved for one interaction point after one year of data taking: 35% at the International Linear 
Collider [^, 8% at the Circular Electron-Positron Collider (sTl and 3% at the Future Circular 


Collider in the electron-positron mode 52 


4 Z 2 -odd Scalar Sector 

4.1 Physical Field Content 

Analogously to the previous section, we define the physical DM fields Aj, Aj and Af, i = 1, ...,3 
through the eigenstates of the corresponding mass matrices of the Z 2 -odd fields rj, fj, rj^ as defined 
in eq. (§. Note that the mass eigenstates are admixtures of both rj and fj. 

We verified that no new gauge vertex is generated through the unitary transformation into the 
physical basis. Furthermore, since the gauge interactions of the r/ and fj are identical up to the 
complex phase i, the total annihilation cross section for two Z 2 -odd into two Z 2 -even fields remains 
unchanged. For r = 0 and with the VEV alignment from eq. (©, the mass eigenvalues for all 
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SU{2)l components of the three Z2-odd scalar fields are given by: 

+ (— Re[(T 3 ^(i)] + 2 cri ( 1 )) v'^ , 

= 2^1 + - (Re[(T3^(i)] + 4 cJi_(i) + 2(T2 ,(i) + 2(T4^(i)) v^. - - 6 m\ , 
= 2^1 + - (Re[(T3^(i)] + 4 cJi_(i) + 2(T2 ,(i) + 2ct4^(i)) v^. + - 6 m\ , 

i 


with 5m\ = v‘^^3 Im[cr 3 ^(i)] ^ + 4 — 0 ' 2 ^(i)cr 4 ^(i) + j. In the limit of r = 0, the eigen¬ 

values from the mass matrices for the Z 2 - and electrically charged scalar fields and the two neutral 
fields are identical and they form complete SU{2)l doublets, i.e. all three components are a gen¬ 
eration of mass-degenerate fields. The couplings to the field break this degeneracy, however. 
The expressions are too cumbersome to be displayed here in a useful way. Fortunately for the 
analysis to follow, it suffices to state that the masses of the fields Ai^ Ai^ Af with i = 1,2,3 are all 
different, and that the relevant mass differences are of the order of uew- 


4.2 Effective Parametrisation 


As can be seen in eq. (43), the masses of the three Ai fields can be chosen independently. The 
same is true for the mass splittings between each generation Ai and its “siblings” Ai, Af, which 
is induced by the (small) VEV of We will thus use the three physical masses and the 

six mass-splittings between the Ai and the Ai, Af, respectively, as free parameters. We define the 
latter for z = 1, 2, 3 by: 


m — m. 


m ; 


— m. 


pi = 


m, 


pi = 


m, 


(44) 


Other important parameters are the effective couplings between the Ai, Ai, Af and the SM-like 
Higgs field Hi. We define this effective coupling for Ac 


d^V 

^ {dAiYidHiY ’ 

with V being the scalar potential from eq. (§. Analogously we can define the couplings w*, 
Especially u>i and are related to the mass splitting parameters pi and pi, respectively, because 
they stem from the same term in the potential but they have no effect on the phenomenological 
study in the following. 

We note the gauge interactions are identical for the Aj=i_ 2,3 and Aj=i_ 2,3 up to a complex phase 
and the couplings between the Aj=i^ 2,3 or the Ai=i^ 2,3 and the Hi depend on the parameters Ki in 
which are independent from the parameters that generate the masses Therefore the 

choice of a specific DM candidate from the three generations is without loss of generality. 

^As mentioned already, the resulting expressions are very cumbersome and not very enlightening, which is why 
we do not to present them here. 
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4.3 Phenomenology 

In this section we will explore the dark matter phenomenology of the above presented model with 
the effective parametrisation from section [4^ In the limit of r = 0 there are two mass-degenerate 
scalar DM candidates in generation i, given by the electrically neutral fields Ai and Ai, both with 
mass and for Depending on the contribution from ueWj two scenarios are given, 

the first one being defined by < m^_ and (scenario 1), and the second one being 

defined by rriA. > mi and < itia , (scenario 2). The fields, the ^i=i 2 3 and the ^*=123 

are linear combinations of the real and imaginary parts of 2 3 ’ with the gauge couplings being 
identical up to a factor i that is associated to the ? 7 i=i, 2 , 3 - The gauge couplings of the ^*= 1 , 2,3 
and ^i=i, 2,3 are thus identical up to a complex phase. In the following we study the limiting case, 
where ^*= 1 , 2,3 consists exclusively in 771 = 1 , 2,3 fields, while ^i=i, 2,3 consists exclusively in 774 = 1 , 2,3 
fields. By a suitable mixing of real and pseudo-real fields a continuous interpolation between the 
two scenarios should be possible. 

We identify the DM field with the electrically neutral Z 2 -odd field with the smallest mass and 
denote it by 4 .dm- The other neutral field will be denoted by in scenario 1, ^dm = Ai while 

in scenario 2, Adm = Ai. Furthermore we introduce the parameters mDM, Pdm,/Odm, which are 
given in the two scenarios by 


scenario I: moM = niAi , POM = 




mAi 


, /5dm = 


- ^A, 


ruAi 


scenario 2: moM = , /5 dm = 


m 


, /5dm = 


"'A. - ^A^ 


(46) 


Ai 


m 


Ai 


Note that in scenario 1, /idm and pdm coincide with the definitions for pi and pi in eq. (44). We note 


here, that even though moM is proportional to a cancellation of terms allows the mass-scale to 
be below the TeV scale despite the assumption in eq. (Q. 

Despite oc ueW) we assume in the following that all the scalar non-DM fields have 

masses above 2mDM- This is done in order to avoid a distortion of the total DM annihilation cross 
section by additional scalar resonances. In the present model and due to the total decay width being 
suppressed by the small lepton-Yukawa couplings, for instance the H 2 resonance would manifest as 
a very narrow peak in fig.[^ leading to smaller minimal values for Ddm at 55 ^dm = 577/72 However, 
since the quark Yukawa interactions are not defined in our model, neither the exact shape nor the 
magnitude of the scalar resonance is properly defined. 

The right-handed Majorana neutrino N[i is also a viable candidate for DM. The combination of 
the Yukawa coupling ip, and the Majorana mass are required to match the light neutrino masses. 


In ref. 16 it was found that for about the Yukawa coupling of the electron, Mjv = O(TeV), 
while for larger values of t/i/, M^r needs to become larger as well. Therefore we shall always keep 
Mtv > 777dm in the following. 


4.3.1 Constraints from Terrestrial Experiments 

The experimental bounds on the masses of charged and neutral scalar fields from LEP-I, which 
entail essentially the width of the Z boson, are given by [^: 

"^A- A A± > 777^/2 (LEP-I). (47) 

Z 5 11 
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Those bounds are stringent on all models that have particles with couplings to the Z boson. 
Furthermore, lower bounds on the masses of charged and neutral scalar particles, denoted by 
and S, respectively, which decay into SM leptons have been derived from LEP-II data [^: 


ms^S± > 80 GeV (LEP-II), 


(48) 


The unbroken Z 2 symmetry that stabilises the DM candidate also forbids the decays of the fields 


Ai, Ai and Af 


into SM leptons, such that the bound in eq. (48) does not apply, 


The OPAL collaboration has published exclusion bounds at 95% conhdence level for supersym¬ 
metric chargino and neutralino particles at LEP-II, see e.g. refs. [53[|54| , at centre of mass energies 
from 189 to 209 GeV. The bound can be expressed as 


cr(e 


+ 0- 


AfAf) < 0.1 pb 


for 


/s = 189 GeV. 


(49) 


Using GalcHEP [40] we get for = mw a cross section for the above process of 0.16 pb. 

The bound in eq. (|49p can therefore be interpreted as a lower bound on the charged DM mass 
m.± > 120 GeV for = ttidm- However, it is obvious from the exclusion plots in 53 that for 
small mass splitting, i.e. poM — 1 the above bound does not apply. In this corner of the parameter 
space the decay products of the charged become very soft and can escape detection due to the 
selection cuts from the experimental analyses. Furthermore, assuming that a number of observable 
events are being produced, the cross section is too small to produce a significant excess over the 
SM backgrounds. Subsequently we will allow for poM and pdm up to 1.2. This leads us to the 
two distinct possibilities for DM with masses below mw- either the mass-splitting between the 
neutral field Adm and the charged held A^,, needs to be large enough such that m .± > mw, or 

UIVI ^DM 

Pdm, Pdm < 1-2. 


As discussed in section 3.2.3 the T parameter is very sensitive to additional scalar doublets. 
Among the “dark sector” helds the mass-splitting parameters pi, pi can be set such that the summed 
contributions of the second and third generation remains compatible with the experimental bounds. 
The mass-splitting parameters pdm and pdm are very important to control the produced relic 


density, as we will show below. We will use eq. (36) therefore to constrain the parameter space 
automatically, unless stated otherwise. 

The constraints from experimental tests on lepton flavour violation also apply to the helds from 
the “dark sector”. In particular the Z 2 -odd scalar helds contribute to the process p —)• ey and are 
thus subject to the strong constraint from the MEG collaboration, cf. hg. (b). Compared to 
the Higgs-like sector those constraints affect the combination of the masses ^±, the Yukawa 

coupling piy and the Majorana mass M^r. In practice given that we are assuming ~ pe and Mm ~ 
TeV , and since these quantities are only constrained through the observed light neutrino masses 
(which gives a bound on the combination), it is not possible to derive meaningful bounds on the 


masses m 




through the bound obtained from the MEG constraint alone. 


4.3.2 Dark Matter Abundance 

The relic density in this section as well as the DM-nucleon-interaction cross sections in the next 
section are computed with the program MicrOMEGAs 56,57 . For all the numerical studies in 


this section we use the Tool for Parallel Processing in Parameter Scans (T3PS) [58]. 
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Figure 5: Numerical results for the DM relic density over the DM mass tudm from a scatter 

scan over the mass and mass-splitting parameters of the Z 2 -odd fields. The considered parameter 
space is defined in eq. (50). The dark green line denotes the observed DM relic density 55 . The 
light and the dark green (non-dashed and dashed) areas denote the scenario 1 and 2, where the 
DM is given by the lightest Ai and Ai, respectively. The constraints from the T parameter on the 
mass-splitting are applied in both cases. In the figure, ^ ^ ttidm is assumed. 
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In order to assess the limits for the two scenarios in our model, we compute the relic density 
rioM for the scenarios 1 and 2 and for the following parameter ranges: 


40 GeV < rriDM < 2 TeV, 

1 < /5dm, Pdm < 1-2 , 

m-DM < 

2mDM < , 

wdm = 0, 

Mat = 10 mDM • 

The resulting values for are shown in fig. by the light green/solid area for scenario 1, 

respectively by the dark green/dashed area for scenario 2. The green horizontal line in the figure 
denotes the presently observed DM relic density [55] . The figure shows that in both scenarios 
parameter values exist, such that the relic density Ddm rnay match the observation for the mass 
ranges below 100 GeV, and between 500 GeV < moM <1-5 TeV. 

We consider the mass ranges for Adm and Adm, where the abundance constraint can be met, 
and conduct two new parameter scans for each scenario with a higher resolution around the areas 
of interest. We set the masses of the i = 2, 3 DM generations to 3 TeV and the DM-Higgs coupling 
wdm to zero, which is identical to uji and Di in scenario 1 and 2, respectively. In fig. by the 
green areas, we show the (moM, Pdm) and (moM, /Odm) parameter space that is compatible with 
the observation of the relic density. The dark/dashed and light green/solid areas denote scenarios 
1 and 2, respectively. 

For the mass mDM ~ 1-4 TeV, the electroweak processes of the form Adm ^dm —^ W~^W~, ZZ 
provide exactly the right cross section to satisfy the abundance constraint. This demands the co¬ 
annihilations to vanish, which in turn implies mass splitting pdm > 1-2. In the figure this results 
in a confined mass range for the DM particle, which is comparable in both scenarios. However, 
since Ddm grows proportional to mDM, the abundance constraint can be met by increasing wdm 
for mDM >1-4 TeV. The DM candidates from the two scenarios considered here have a mass that 
is confined to one of the two intervals: 


m^i G 


47 GeV to 74 GeV 
600 GeV to 3.6 TeV 




50 GeV to 74 GeV 
620 GeV to 3.6 TeV, 


(51) 


where the upper bound on the mass comes from the perturbative limit LOi, uit < 1 for all i, and 
requires the condition pDM > 1-2 and /5 dm > 1 -2, i.e. suppressed contributions from the other “dark 
sector” fields. 


4.3.3 Experimental Limits from Direct Dark Matter Searches 

For /Odm,/5dm > 1-2, the observed value for the relic density can be used to constrain the ^dm- 
Higgs coupling, parametrised by wdm- We analyse both scenarios for the two values of pdm = 2.0 
and /Odm = 1-15 and we set the mass splitting between the two neutral components /5 dm = 2.0. 
All the other masses are set to 10 x mDM, to avoid co-annihilations. We remark that regarding 
the numerical computation of the relic density both values of /?dm = 2 and /Odm = 10 lead to the 
same result. This is due to the Boltzmann suppression of the “dark sector” fields, which effectively 
removes their contribution to the total annihilation cross section for mass-splittings much larger 
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Figure 6: Left: The relative mass-splitting parameter pdm (as defined in eq. (46)) over the DM 
mass. The light and dark green (non-dashed and dashed) areas denote resulting values for the 
relic density Ddm in agreement with observation for scenario 1 and 2, respectively. Right: Mass 
splitting parameter pdm between the two neutral DM fields Ai and Ai for niAi < over the 

DM mass. The colour-code is the same in both panels. The constraints from the T parameter on 
the mass-splitting are applied in both cases. In this figure, mHi ^2 > ^^i >2 ~ ^ been used. 




(a) scenario 1 (b) scenario 2 


Figure 7: The left panel denotes scenario 1, where Adm = Ai, tudm = ^Ai and wdm = while the 
right panel denotes scenario 2, where ^dm = Ai, tudm = a;DM = for ruAi < rriAj 




. In 


both panels the green and grey bands denote the mass-splitting between Adm and = 2.0 

and pdm = 1-15, respectively. Parameter values inside the green and grey bands result in a relic 
density Ddm that matches the observation. The red area denotes the exclusion constraints from 
the direct DM search experiment LUX 


59 . In this figure, p = 2.0 and all the other masses are 


chosen to be above 100 GeV. 
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than 20%. The particular setting of the masses is thus somewhat arbitrary since it has no further 
impact on the result. Note that the “corner” of the parameter space considered here does not lead 


to tension with the bounds on the T parameter in eq. (34). 


The resulting constraints for wdm are shown in the two panels in fig. where the green and 
grey band denote pdm = 2.0 and pdm = 1-15, respectively. We observe that the allowed region 
for cudm is being “squeezed” to smaller values, for poM approaching the critical value, which is 
given by the lower bound in fig. Also displayed in the figure are the LUX constraints, which are 
denoted by the red area, which affect the spin-independent ADM-nucleon cross section, The 

experimental upper bound is of order zb (10® zb = 1 fb or zb = 1000 sheds) for tudm between 10 
and 100 GeV, and order 10 zb for tudm ~ 1 TeV. It is interesting that fig. (b) shows the 
possibility to generate the correct DM relic density for scenario 2 also for ttidm > j which is due 
to a destructive interference in the total annihilation cross section that is not present in scenario 1. 
Note, that the LUX bound relaxes for larger masses, such that for ruyi. beyond 3 TeV a coupling 
|<^dm| ~ 1 is possible. 

The fact that the relic abundance can be generated in our model via the co-annihilations is 
visible in the “squeezing” of the contour for wdM; which occurs when pdm becomes smaller than 
1.2. Therefore the direct search constraints yield strong bounds on the DM interaction with the 
Higgs boson, but do not exclude the possibility of DM with 47 GeV < tudm < 74 GeV. 

The spin-dependent interaction between the DM and the nucleon, with the corresponding cross 
section is mediated by the weak gauge bosons. Experimentally, (7^^ is constrained especially 
by IceGube (60] w ith upper limits of 0.1-1.0 fb, followed by the upper bound of 4.5 fb from Super- 
Kamiokande [6l|, both valid for the mass range from ~ 100 GeV to 1 TeV. For masses below 100 
GeV, the best upper bounds to the spin-dependent DM-nucleon interaction is given with ~ 4 fb 
from COUPP 62 and SIMPLE 63 . The interaction between the field Adm and the nucleon takes 


place only at the loop level via exchange of at least two W ov Z bosons, and the cross section is 
suppressed by compared to the usual tree-level expression involving the exchange of a single 
vector boson (^, which yields a typical cross section 0(0.1) fb, for weak cross sections 0(pb). 


4.3.4 Signatures for Dark Matter at Hadron Colliders 

Direct searches for signatures of dark matter with charged components at the LHC have been 
considered e.g. in refs. |65l|66|, and very recently also including the hadron-hadron mode of the 


Future Circular Collider (FCC-hh) planned at CERN 67 


The search for a signal in LHC data is very promising, provided the DM mass is within reach. In 
particular, event topologies with an imbalance in transverse momentum and one (or two) hadronic 
jet(s) or photon(s), respectively, may show an excess over the SM background. In fact, for an 
appropriate choice of masses it might be possible to explain observed excesses in events with 
opposite sign di-lepton, jets and missing transverse momentum by ATLAS and CMS - 
due to cascade decays of the three generations of Z 2 -odd scalars. Furthermore, due to the specific 
vacuum alignment of our model, the interactions between fields of the second generation with those 
of the first and/or third generation are suppressed by the small parameter r, such that the second 
generation fields can have a long life time. This in turn allows for signals with displaced vertices, 
which benefit from a very low SM background, which was recently considered in ref. 67 in the 
context of an effective field theory approach. 

Recently, the CMS collaboration has also published bounds on the spin-dependent-DM-nucleon 
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interaction cross section from an analysis of 19.7 fb ^ LHC data at ^/s = 8 TeV 70 . Therein the 


collaboration considered on the one hand an effective operator approach, which is not applicable 
to our model. A second part of the analysis considers a s-channel mediator, which couples in a 
specific way with quarks and the DM particles. In considering the s-channel mediator being given 
by the Z-boson, constraints on the present model parameters can be obtained from the LHC data 
in a similar analysis. Such an analysis is beyond the scope of the present paper. 


5 Conclusions 

In this paper we considered an A 4 x Z 2 x Z^-symmetric extension of the Standard Model (SM), 
wherein the S'[/(2)i-doublet leptons are embedded into A 4 -triplet helds and the matter content 
also includes a Z 2 -odd Majorana neutrino. The scalar sector consists of two A 4 -triplet SU{2)l- 
doublet fields $(^=^’ 2 ) among which one is charged under the Z 2 symmetry, and further a SM-singlet 
A 4 -triplet Z^-odd scalar field, whose vacuum expectation value (VEV) breaks A 4 and Z^. This 
renormalisable model has previously been shown to successfully account for the charged lepton and 
neutrino masses and mixings through a specihc alignment of the VEVs. The small masses of the 
active neutrinos are generated with a radiative seesaw mechanism. The scale of the A 4 and Z 2 
breaking VEV can be as small as O(TeV) and thus remain comparatively close to the electroweak 
scale. We presented the A 4 x Z 2 x Z^-invariant scalar potential for the considered field content 
and discussed under which circumstances the specihc VEV alignment - which gives rise to large 
leptonic mixings - can be realised in this model. 

An interesting feature of the model is that the second and third generation Z 2 -even scalars 
have couplings to leptons that are exclusively lepton havour violating (while the hrst generation 
of the Higgs-like, Z 2 -even scalar helds, has coupling strength to leptons that is very close to the 
corresponding coupling of the SM-Higgs boson). We have discussed phenomenological constraints 
on the masses and mass-splittings for the extra Z 2 -even scalar helds from the oblique parameters 
and direct searches for lepton havour violation at low energy and at LEP. The strongest lower bound 
on the mass scale is 140 GeV due to the constraints from the MEG collaboration on /x —)• ey. At the 
LHG, this model should have created around 0(10) lepton havour violating events in the current 
data such as /urjj and erjj, where j is an energetic jet, and less than 0(1) events of efijj events. 
It is possible to interpret the excess in /rr events found by the CMS collaboration as being caused 
by the helds H 2 , H 2 , which implies that the MEG collaboration should observe /r —)■ ey very soon. 
Alternatively, a perturbation of 0(0.1) of the VEV alignment can lead to a LEV branching ratio 
of the SM-like Higgs boson Hi —)• /rr that can explain the observed excess at the LHC. 

The Z 2 -odd held content of the model yields a “dark sector” with the DM candidates being 
the Majorana neutrino and the neutral components of the A 4 -triplet S'f7(2)i-doublet among 
which we investigated the properties of the scalar DM candidate. First we discussed the existing 
limits from LEP I and H. In particular, the OPAL collaboration provides constraints on the masses 
and mass splittings for scalar DM with masses below 100 GeV, which turn out to be negligible when 
considering the DM properties of the model. Then we derived the constraints on the model param¬ 
eter space from the observed dark matter relic density and from direct search experiments. The 
experimental constraints from the LUX collaboration imply that the DM particle barely interacts 
with the Higgs boson when its mass is 0(100) GeV. 

We found that it is possible to account for the observed relic density for DM with a mass in the 
interval between 47 and 74 GeV or in the interval 600 GeV and 3.6 TeV. An important parameter 


21 



in the model is the mass-splitting between the neutral DM particle and its charged “sibling”, and 
also the mass-splitting between the three generations of DM particles. 

It is promising to establish the connection between Dark Matter and the scalars at the LHC and 
especially at future lepton colliders by studying the decay properties of the 125.7 GeV SM-Higgs- 
like scalar with higher precision. Last but not least, direct searches in the high-energy data from the 
next run of the LHC and at future hadron colliders might reveal the extended scalar field content 
of this model by detecting di-lepton excesses, or multiple (scalar) resonances for centre-of-mass 
energies of several TeV. 
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A The product rules for A 4 

The following product rules for the A 4 group were used in the construction of the Lagrangian: 

3®3 = 3,©3a©l©l'© 1", (52) 


l©l = l, T©l" = l, T©T = l", T©T = T, (53) 


Denoting (xi,yi,zi) and {x 2 ,y 2 ,Z 2 ) as the basis vectors for two A 4 -triplets 3, one finds; 


(3 © 3)1 
(3 © 3 ) 1 / 
(3©3h// 

(3 © 3 ) 3 ^ 
(3 © 3 ) 3 ^ 


xiyi + X2y2 + xsys, 

xiyi +ujX 2 y 2 

Xiyi +U}‘^X2y2 + WX3?/3, 

(2:22/3 + 2:32/2,2:32/1 + 2:12/3,2:12/2 + 2:22/1), 
(2:22/3 - 2:32/2,2:32/1 - 2:12/3,2:12/2 - 2:22/1), 


(54) 


where a; = The representation 1 is trivial, while the non-trivial V and V are complex 

conjugate to each other. Comprehensive reviews of discrete symmetries in particle physics can be 
found in refs. 
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B Scalar potentials with A 4 triplets 

The most general potential for two fields, charged under separate Z 2 and simultaneously 

A 4 triplets and SU{2) doublets has a large number of possible contractions: 


Vu = Y. 


i=l ,2 L 


^2 + Ai4i) ($W^.«)i(|.«$W)i + A2,(i) ($«$«)i($«4>«)] 
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I-1 
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F- 1 




I-1 


I-1 


I-F 


n I 


I-^ 


F I 


I-F 


F I 


+ Kg + |kio| e‘^i° 

I— I I —I 


I-1 


I-1 


I- 


I-1 


I-1 


I 


+ |ki 3| e'^13 ($(l)$(2))^(|,(l)$(2))^ ^ gi/3i4 ($(1)$(1))^(${2)${2))^ 

I-^ 


F I 


I-1 


I-1 


+ |ki 5| ($(1)$(2))^((J)(2)$(1))^ |^^g| gi/3l6 ($(1)(I)(2))^,(${1)(J){2))^,, 

I-1 I-1 

+ |ki7 | + h.c., 


where is the charge conjugate of lines above helds signify SU{2) contractions 

using the tensor e = \a 2 and subscripts of brackets denote the type of A 4 contraction used for the 
helds within. All couplings are either real or their phase is stated explicitly. All additional SU{2) 
and A 4 invariant terms can be expressed as linear combinations of the above. Assuming manifest 
CP symmetry [74| , all couplings must have vanishing imaginary parts. 

Likewise if we introduce a real scalar held y, which is a singlet under the SM gauge group and 
a triplet under A 4 , we hnd 


^x<t> — X/ 

i=l ,2 L 


I-1 




+ k3,«l +^4,« (^>^'^X)1'(‘J>^*^X)] 


I-1 


(55) 


+ h.c. 


for the potential mixing both sectors and 

^xx =/^x (xx)3,(xx)3s + d 2 (xx)i'(xx)i» (56) 

for the potential only depending on y. In order to make the vacuum as shown in eq. Q an extremal 
point of the potential, the couplings (73 ( 2 ) and 174 ( 2 ) must be absent or negligibly smalj^The mass 
terms for and y fulhl the relations 

hi = “(-^1,(2) + ^2,(2) + "^3,(2))^^^ “ cri,(2)'yx ’ 

®For details on the effects when this is not the case, see Appendix [3.2.2 
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^ ((4di + d2)vl. + 0 - 1 ,( 2 )V^) . 


(58) 




2 

X 


C Solution to the pathological vacuum 


In the sub-potential while all derivatives at the position x = ^le set to zero 

with an appropriate condition for or the mass square eigenvalues for the x subsector are 
then not positive definite. 

This is an issue for the model being considered as it relies on the VEV x = 1), namely 

for obtaining eq. Q. The problem can be solved for instance by allowing r ~ 1 or by a redefinition 
of the VEV direction in A 4 space, to either (1,1,1), a permutation of (1,0,0) or even a complex 
direction (1,0, i). In those cases however, we lose the connection to the flavour model in |^, and 
it remains to be shown that large leptonic mixing angles are still possible. Due to the vanishing 
entry, the latter VEV direction (1,0, i) has comparable phenomenology to (1,0,—1) in terms of 
leptonic mixing angles (cf. eq. Q), but a detailed analysis of this is beyond the scope of the present 
paper. We note that the change of the VEV direction in A 4 space does not affect the DM analysis 
in section 4.3, due to the effective parametrisation used. However, changing the VEV has a strong 
effect on everything else. 

The solution we prefer to use here is to introduce another complex singlet scalar field ^ in the 
1" representation of A 4 . Its interactions to x lead to added terms in the scalar potential 


~ 'I'l ^(xx)i' + 72 ??? + + A^ ^?(xx)i'' + h-c. • (59) 

Together with the terms in l^X’ shown that this potential can lead to (depending on 

the parameters) a VEV {^) = and the aforementioned (1,0,—1) VEV direction for X) with 
positive mass squares of similar magnitude for all mass eigenstates around this vacuum. Since the 
terms ^ and (xx)i" have the same assignment under all symmetries, the presence of ^ also leads 
to additional couplings of form similar to the CT^^(i) term of 1 ^ 0 . However, due to the form of (x) 
one finds that 

((xx)i") = -^ (0 , ((xx)i') = -% , (60) 

which means that these additional couplings between and ^ can simply be absorbed into cr 3 _(j) 
as far as masses for the and the VEV alignment (with potential pertubations) are concerned. 
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